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An approximate general solution of the differential equation 

d 2 w/dz 2 =f(z)ic 

is 

w^Af-" 4 exp (ff l l 2 dz) + Bf~ i l 4 exp (-ffl^dz). 

To represent a particular solution throughout a given complex domain, differing pairs of values of the 
arbitrary constants A and B are generally needed in various subdomains. Strict error bounds are 
established in this paper for the coefficients in the linear equations connecting the pairs of values of 
A and B, in the case in which the z-domain under consideration is unbounded and contains a simple 
zero and no singularities of /(z). 

The results are used to place the "WKBJ" method on a firm mathematical foundation. 

Key Words: Differential equations, asymptotics, phase-integral, WKBJ method, turning point, 
connection formulas, error hounds. 

1. Introduction 

A form of solution of the differential equation 

cPw ri , 

— =f(z)w, (1.01) 

which is frequently used in applied mathematical problems, is given by 

w - 4/*" 1/4 exp ( Sf m dz) + Bf- 1 ' 4 exp (- J f'hlz) . (1.02) 

in which A, B are arbitrary constants. This is commonly referred to as the Liouville-Green approxi- 
mation, the WKB or WKBJ approximation, or the phase-integral solution. A necessary condition 
for the applicability of the approximation is that/~ 1/4 be a bounded, slowly varying function of z [l]. 1 
This condition automatically excludes from the region of validity most kinds of singularity of/, 
and also zeros of/ the so-called turning or transition points of the differential equation. 

A further condition concerns the monotonicity of the function Re ff ll2 dz (see [1], sec. 8, and 
sec. 4, below). This requirement has the effect of dividing the z-region under consideration, D, 
say, into a number of subregions, in each of which a different combination of values of A and B 
in (1.02) has to be used in order to represent the same solution of (1.01). This is, of course, an 
example of the Stokes phenomenon. The linear equations linking the values of A and B in one 
subregion with those in another are called the connection formulas or Stokes equations. The 
purpose of this paper is to determine the connection formulas in the case in which D contains 
no singularity of /(z), and just one simple turning point, z , say. 



1 Figures in brackets indicate the literature references at the end of this paper. 
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Two principal methods are available for attacking this problem. The first, applicable when 
f(z) is real on the real axis, starts with the Liouville-Green approximation in real exponential form 
on one side of z , then traces its behavior along a semicircle in the complex plane to the real axis 
on the other side of zq. The condition that the solution must be real whenz is real suffices to deter- 
mine the main features of the appropriate connection formula. This method was originated by 
Zwaan [2], and has been developed further by Kemble [3] and Furry [4]. 

The second method consists of using approximations to the solutions in terms of Airy functions 
(or equivalently Bessel functions of order one-third), the region of validity of which includes z () . 
Away from Zq the Airy functions are approximated by their asymptotic forms for large argument, 
leading immediately to the Liouville-Green forms. In the original version [5]-[7], the Airy-function 
representations were constructed by replacing /(z) by a linear approximation in the neighborhood 
of z ; this is, perhaps, the true "WKBJ method." In more recent work [8]-[13] more powerful 
uniform approximations in terms of Airy functions have been used to derive the connection formulas. 

Hitherto, the method of Zwaan has not been put on a sound mathematical footing: writing in 
1934 Langer [14] commented ". . . an estimation of the relative error involved in the result is 
almost impossible . . .". The treatment of the Airy-function procedure has been more satis- 
factory. By supposing that /(z) contains as multiplicative factor (the square of) a large parameter 
a, say, it is shown in the references cited that the connection formulas are asymptotically correct 
as \u\ -» °°. 

In the present paper we give rigorous analytical derivations of both methods, and supply 
explicit bounds for the errors in the connection formulas (thus overcoming one of Langer's objec- 
tions to Zwaan's method). Since we are dealing with bounds there is no need for the explicit 
introduction of a parameter u in the theory: asymptotic properties in applications can be regarded 
as immediate consequences of the form of error bound. 2 We also dispense with the requirement, 
in Zwaan's method, that f(z) be real on the axis. 

The starting points of our investigation are the theories of error bounds for the Liouville-Green 
and uniform Airy approximations which have been developed by the writer in [1] and [15]. 

Since the present paper and its sequel [19] were written, there has appeared a monograph on 
the same subject by N. Froman and P. O. Froman [20]. Although the objective of these authors 
is essentially the same as that of the present writer, the methods used and results achieved are 
sufficiently different to warrant publication of this paper and [19] in substantially their original 
forms. Some comparisons with the results of [20] are given in section 11. 

2. First Error Theorem 

We begin with the differential equation 

^ = {1+/U)}w, (2.01) 

in which j{z) is a regular (holomorphic) function in a simply-connected complex z-domain D. 

THEOREM 1. Equation (2.01) has solutions Wj(z) 0=1? 2) regular in D, and depending on 
arbitrary fixed points aj o/D, such that 

w j (z)=e ±z {l + e j (z)}, (2.02) 

where 

\ej(z) | , ]e'j (z) | , \ej(z) ± ej (z) \ ^ exp {V #j (C) } - 1 (2.03) 

whenz E H,-. 



2 In the course of a valuable historical introduction to phase integral methods, Heading ([13], p. 8) critizes other writers, particularly Budden [12], for failing to 
standardize a large parameter in their analysis. The theory of error bounds presented in this paper provides a broader basis for justifying the work of these writers. 
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In this result, the upper signs are to be taken when j = 1, and the lower signs when 7 = 2. The 
regions of validity H, (/==1, 2) comprise the aggregate of those points z for which there exists a 
path^j joining z with a, and lying entirely in D, with the properties: 

(i) 0>j consists of a finite number of Jordan arcs, each having a parametric equation of the form 
t = t(r) with t"{r) continuous and t'(r) nonvanishing. 

(ii) Re t is monotonic nonincreasing (/= 1) or monotonic nondecreasing (/=2), as / traverses 
&j from z to dj. 

Lastly, f]^ (G) denotes the variation along 3L$ of the function 

G(z)=\ f(z)dz; (2.04) 

thus 

^(G)=J \f(t)dt\. (2.05) 

Theorem 1 is essentially a special case of Theorem 7 of [1]. We record the main steps of the 
proof, however, since they are needed in the derivation of further properties of the error terms in 
the next section. 

For illustration, we takey= 1. From (2.01) and (2.02) we find that 

e;'u)+2€;u)=/u){i+€,u)}. (2.06) 

Therefore €1(2) satisfies the integral equation 

6,(2) - •> Ml- c*'"*>} { 1 +€,(/)}/( t)dt, (2.07) 

in which the path of integration is taken to be & \. The solution of this equation is constructed 
as a convergent series 

€1(2) = 2 /i,„UK (2.08) 

in which l\ % 0(2)= 1, and 

/.. „(z)=l \ { l-e*'-*>}j\t)l UH -i(t)dt (n ^ 1). 

J (h 



(2.09) 



In consequence of Condition (ii) on & \ we have 

|e 2( '-*'| « 1, (2.10) 

and by substituting this bound in (2.09) we readily establish by induction that 

\l,.n(z)\^{i / ai ,AG)}"lnl (»>0), (2-11) 

showing that the series (2.08) is uniformly convergent. In order that solutions of the integral 
equation (2.07) also represent solutions of the differential equation (2.06), some restrictions are 
needed on the path &\. Those given in (i) above suffice (compare [16], pp. 491-492). 

The first of the inequalities (2.03) now follows by summation of (2.11). The other two are 
derived similarly from the differentiated form of (2.09). 
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3. Further Properties of the Error Terms 

For conciseness, we shall refer to a path &j which fulfills Condition (i) of section 2 and is such 
that Re t is monotonic as t traverses 0> - } as a progressive path. If the monotonicity is strict every- 
where, we say that the path is strictly progressive. It may be noted, incidentally, that Condition 
(i) is quite weak, and can be fulfilled, for example, by any finite chain of straight lines and circular 
arcs. 

The condition that D be simply connected, mentioned in the opening paragraph of section 2, 
merely ensures that each Wj(z) is single-valued in D. For a multiply connected domain, each 
branch of Wj(z) satisfies (2.02) and (2.03) whenever a progressive path can be found in D joining 
z with a,j. 

In section 2 we tacitly assumed that the points aj are finite. Theorem 1 still holds, however, 
if D is unbounded and the cij are points at infinity on progressive pathsi?;, provided that 0* j coin- 
cides with J2j as cij is approached and that the variations (2.05) converge. The convergence of the 
infinite integrals (2.09) is then an immediate consequence of the inequality (2.11). 

A situation of importance in the later sections occurs when d\ and a 2 are distinct points at 
infinity, and can be joined by a progressive path > lying in D. From (2.03) it is seen that 6i(z) and 
e[(z) both vanish as z— > a\ along &\ indeed we may regard these properties as the boundary con- 
ditions defining W\{z). What happens, however, to €\{z) and e[(z) as z— » ap. 

THEOREM 2. //]Re z|— > °o as z— > a\ and as z~> a?, then 

€i(z) — » a constant e\ (a 2 ) , say, and e[(z) — » as z—> « 2 ,] 

(3.01) 
£2(2) "~ > a constant ei(ai), say. and €9(2:)— » as z—> a u ) 

ei(a 2 )=€ 2 (ai), (3.02) 

and 

|€i(<b)| ^Hexp {/>,(G)}-1|. (3.03) 

To prove these results, let 

a u „(z)= I f(t)lu„-i{t)dt, Pu„(z)= J e 2it - z) f(t)l Ull -i(t)dt. 

J<i] J(i\ (o.04) 

Then from (2.08) and (2.09), we have 

6i(z)=i la Ull (z)- h Sj8i.„(z), e((z) = 20i,„(z), (3.05) 



con- 



the summations extending in each case from n — 1 to rc = °°. Since by hypothesis I \f{t)dt\ 

J a 1 

verges, it follows that for any given positive number 8 which is less than the value of this doubly 
infinite integral there exists a point z on such that 



\ <l2 \f(t)dt\ = 



8. (3.06) 



Let z be any point of between z and a 2 (fig. 1). Then by splitting the integration range in the 
second of (3.04) at z', and applying (2.10), (2.11), we find that 

n\ (n-l)l 



\Pi.M\ * |«**-.>| lila^ai + U < :: «w - 8 . {301) 
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Figure 1. 



Therefore 



\XPu„{z)\ ^ \e 2(z '- z) \ [exp {/ ril ,*<G)}-l]+8exp {/ flI ,„ t (G)}. 



(3.08) 



With the assumed conditions, e 2(z '~ z) vanishes as z-> a 2 . Hence 2/8i,»(z) also vanishes, that is, 
e;(z)->0. 

Next, we have 



\a\, „iz) — a], ,\z 



/; 



f{t)lun \(t)dt 



1 « U AGJ\»-{1 ( , U AG)}» 



(3.09) 



ISai.wUHSai.Hte')! ^exp { ■/ „„*(G)}-exp { i W| ,*KG)}. 



(3.10) 



The righl side of tliis inequality vanishes as z and z' lend independently to a->* therefore Sai, w (z) 
tends to a constant limiting value. The proof of the first line of (3.01) is now complete. The 
second line follows by symmetry. 

To prove (3.02), we construct the Wronskian from (2.02): 



^(w 2 ,w 1 )= {1 + €2(2)} {l + €i(z) + e,'(z)} + {l + €1(2)} {l-r-e 2 (2)- e! 2 (z)}. 



(3.11) 



From (3.01) we see that this reduces to 2{l-f €2(^1)} at z — au and 2{1 -f €1(^/2)} at z= a 2 . Since 
the Wronskian is constant. (3.02) follows. 

Lastly, we derive from (2.11) and (3.04) 



\a 1 ,n(a 2 )\^{r ait a 2 (G)} n lnl 



3.12) 



The improved form (3.03) of the main result (2.03) for the special value z = a-i now follows imme- 
diately from (3.05). 



4. Basic Error Theorem for the Liouville-Green Approximation 

We now proceed to apply the results of the preceding two sections to the differential equation 
(1.01). A useful generalization is achieved on replacing the coefficient of w in this equation by the 
sum of two functions /(z) and //(z), say. In applications the second function h(z) will either be 
zero, or small, in some sense, compared with fiz) (cf. [1|. sec. 4(h)). From the standpoint of the 
error analysis, however, the relative magnitudes of f(z) and h(z) are immaterial. 

Thus we consider the differential equation 



d 2 w 
dz 7 



{f(z) + h(z)}w, 



(4.01) 
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in which f(z) and h(z) are regular in an unbounded domain D, not necessarily simply connected, 
and f(z) is nonvanishing. The Liouville transformation 



(4.02) 



£ = £(*)= f{f(z)yi*dz, r={/(z)}^V 
throws (4.01) into the form 

d 2 W/de= {l + <p(QW, (4.03) 

in which 

, n _h(z) Wz)f'(z)-5{f'(z)} 2 h 1 rf» / 1 \ (404) 

Equation (4.03) is of the form (2.01). The z-domain D is mapped on a ^-domain A, say, and (p(g) 
is regular in A. By applying Theorem 1 and then restoring the original variables we establish the 
following: 

THEOREM 3. Let z — c be an arbitrary point at infinity in D such that |Re f(c)| =°°, and let 
the branch of i;(z) be chosen in such a way that Re g(z) is positive as c is approached. Then the 
differential equation (4.01) has a regular solution, given by 

w(z) = {f(z) }~^e-^{l + e(z) }, (4.05) 

where 

|e(z)|, |/-»'V(z)|, |e(*)-{e' (*)/£' &)}|<exp {f..e(*)}-l (zGH(c)), (4.06) 

F<d=f V (£)di^f\p~£ {$}*. (4.07) 

The region of validity H(c) comprises those points z which can be joined in D to c by a "£ -pro- 
gressive" path ^, that is, one whose map in the £-plane is progressive (sec. 3). All branches in 
(4.05) to (4.07) are continuous along £? , that of / 1/2 being the square of/ 1 / 4 . The variation in (4.06) 
is taken over $P\ naturally we assume that it converges at c, otherwise (4.06) is meaningless. 



5. Fundamental Connection Formula for a Simple Turning Point 

In this section we suppose that/(z) and Hz) are regular in an unbounded, simply connected 
domain D, and that/(z) has just one zero in D, a simple zero at an interior point zo, say. We 
define tj(z) by (4.02), fixing the arbitrary lower integration limit at z : 

£»= P if(t)} ll2 dt. (5.01) 

Clearly g(z) is an analytic function in D whose only singularity is a branch point at z . If the 
Taylor expansion of f(z) in the neighborhood of z is denoted by 

f(z)=^f n (z-zo) n , (5.02) 

71=1 
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then 



t{z) = -f\Hz-z () fl*{\+0(z- 



■*>)}. 



(5.03) 



The set of points in the z-plane for which 



Re f(z) = 



(5.04) 



we call the principal curves 3 associated with the turning point z - Using the theory of conformal 
mapping, we see that each principal curve is a regular Jordan arc which can terminate only at z 
or at the boundary of D. Moreover, no principal curve can intersect itself or any other principal 
curve except at z<). Equation (5.03) shows that three principal curves emerge from zo, inclined 
there at angle 2-7r/3 with each other; we denote them by #1,^2,^3, enumerated in the positive 
rotational sense. These and other principal curves divide D into a number of subdomains. Those 
with 2o on their boundary we call the principal subdomains, and denote them by Di, D2, D3, as 
indicated on figure 2. Clearly each D, is itself simply connected. For terminological convenience 
in what follows we shall enumerate j with modulo 3; thus D = D.j, D 4 = Di. 

We define £j (z) (/= 1> 2, 3) to be that branch of £ (z) which is continuous in D cut along the cor- 
responding^, and whose real part is positive in D, and negative in the other two principal sub- 
domains. Analytic continuation across the cut^j is expressed by 



6{(z-zo)e ±2 -}=-6(z-z ); 

compare (5.03). We define g(z) to be the four-valued function 

giz)={f(z)}->l\ 
Its analytic continuation is expressed by 

g{(z-z {} )e ±27Ti }= Tig{z-z {) ). 



(5.05) 



(5.06) 



(5.07) 



We denote by g/z) (/=1, 2, 3) an arbitrary constant multiple of any branch of g{z) which is con- 
tinuous in D cut along the corresponding^,. In order to achieve a symmetric form of final result, 
we stipulate that these multiples must be chosen in such a way that 



gj (z) = e™l« gj -Az) (z G D,_, U Dj) 



(5.08) 



ed anti-Stokes lines [12], [13]. 



Figure 2. Principal I subdomains Dj. 
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The self-consistency of this equation for j= 1, 2, 3 follows from (5.07). 

Let us suppose that each D, is unbounded, and that Cj is an arbitrary point at infinity in D, 
such that | Re f(c/)| = °°. Excluding the point zo from D, and then applying Theorem 3 of section 4, 
we find that eq (4.01) has solutions wj(z) (j= 1, 2, 3), given by 

Wj(z) = gjiz)e-¥ Z) {I + ej(2)}, (5.09) 

where 

|€j(z)|J/- 1 / 2 €;U)|j€j(2:)-/- 1 / 2 €;(z)| ^exp {i z . r .(F)}-l (z E H(c/)) (5.10) 

anql F is defined by (4.07). In (5.10) the branch of/~ 1/2 is l/£j(z), and H(cj) is the set of z-points 
which can be joined to Cj by a ^-progressive path &j lying in D. (We observe, in passing, that the 
corresponding principal curve ^j is excluded from H(c/), since passage of .^/ through z is not 
permitted.) *f z, Cj (F) denotes the variation over &j of F. 

Although each solution Wj(z) satisfies (5.09) and (5.10) in only part of D, the theory of linear 
differential equations shows that each Wj{z) exists and is regular throughout the whole of D. Fur- 
ther, there is an identity of the form 

AiWi(z) -\-A 2 w 2 (z) +^3^3(2) = 0, (5.11) 



in which the Aj are independent of z. Equation (5.11) is the connection formula for the turning 
point Zo, and the values of the Aj constitute our primary quest. In order to obtain them, we assume 
that ^-progressive paths exist in D linking each pair of reference points c- h c/,. 

Consider (5.11) in Dj. Substituting by means of (5.09) and (5.08), and using the relations 
£2(2) = f 3(2) = — £iU), we derive 

A 1 e-fM{l + e 1 b)}+A 2 <pM** g Hl+e s b)}+ (5.12) 



Now let z— » C\. Then €\{z) vanishes and, by Theorem 2, e 2 (z) and e-Az) tend to constant values e-zici) 
and e^ici). Also by hypothesis Re fi(z)— » + °c. Hence we have 

e 7n76 {l + e 2 (ci)}^2 + e- 7rj '/ 6 {l-r-€3(ci)}^ 3 = 0. (5.13) 

Similarly, 

e -mi6{i + €l (cd}Ai + e^Hl + e 3 (c 2 )}^3=0, (5.14) 

e 7n/6 {l + e,(c3)}^i + e-^ 6 {l + e 2 (c3)}^2 = 0. (5.15) 

The consistency of this set of homogeneous equations for the Aj immediately follows from the 
relation Cj(ca-) = €a-(cj) obtained from eq (3.02). Solving for the Aj and substituting back in (5.11), 
we arrive at the desired form of the connection formula: 

{l + e 2 (cs)}^Mz) + e- 2m7 Hl + ^ (5.16) 
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The error terms are hounded by (compare (3.03)) 

|€/a.)|«i[exp{'/ ,,,,(F)}-1], 



that is, 



\6jiC k ) 



■ exp 



pit' 



_L^1 

V ,/4 dz 2 



\dz\ 



(5.17) 



in which the variation (or integral) is evaluated along a f -progressive path linking Cj with o,. 

We summarize this section so far. The assumed conditions on f(z) and h(z) are stated in the 
opening sentence. The principal subdomains D y associated with the turning point z () are defined 
in the second paragraph. Three solutions wj(z) (/=1, 2, 3) of the differential equation (4.01) are 
then constructed, given by (5.09) and (5.10). Here £,(z) is the branch of the integral (5.01) which 
has positive real part in D y ; gj(z) (j= 1, 2, 3) is any set of constant multiples of branches of {/(z)} _1/4 
satisfying (5.08); Cj is any point at infinity in D, such that |Rc £(c/)| = °°. The error bounds (5.10) 
are valid at all points z which can be joined to Cj by a ^-progressive path lying in D. Finally, pro- 
vided that the cj can be joined in pairs by paths of this kind, the fundamental connection formula 
between the Wj{z) is expressed by (5.16) and (5.17). 

Remark. It will be observed that an outstanding problem in the application of the connection 
formula (5.16) will be the determination of suitable ^-progressive paths linking together the points 
at infinity cj. This appears to be an unavoidable feature of this method, 4 and its solution will 
generally necessitate a study of the actual conformal transformation between the £ and z-planes. 

Usually it will not be difficult to demonstrate the existence of suitable ^-progressive paths. 
These paths are not unique, however, and the most effective use of the connection formula clearly 
calls for the determination of those paths which minimize the total variation of the function F along 
their length. Roughly speaking, this means that we need to keep the ^-progressive paths as far 
away as possible from all the singularities of ' f(z) and h{z) and zeros of f(z) (including Zo), subject 
to fulfillment of the monotonicity condition. 

The construction off-progressive paths in the neighborhood of infinity is discussed in section 7 
below, and an example of the determination of ^-progressive paths which minimize the variation 
of an analytic function along their length is given in [17], section 7. 

6. Airy's Equation 

The simplest example to which the theory of section 5 is applicable is provided by the dif- 
ferential equation 

d 2 wldz 2 = zw. (6.01) 

In the notation of sections 4 and 5, we take f(z) = z and h{z) = 0. Then Zo = and 



f(z) = |z* 2 , g(z) = z-V\ 



(6.02) 



D comprises the whole z-plane, and the principal subdomains D/ are the sectors of angle 27r/3 
indicated in figure 3. 

We define g\(z) to be the principal branch of ' z 1/4 . Then from (5.08) 

g-Az) = e 7Tilli z- ] l 4 (z ED,U D,), g s (z) = e-™ l6 z- 114 (z G D, U D 3 ). (6.03) 

* Compare |20|. chapter 4. 
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FIGURE 3. Principal subdomains for Airys equation. 
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Figure 4. %-plane. 



We take Ci= + o°, C2 = °° e 27r//3 , C3 = °° e~ 27ril3 . By considering the £-map we see that each region 
H(cj) comprises the part of the z-plane remaining after the removal of the corresponding principal 
curve ^ j. From (4.07) we have 



Fn= 5 = 5 
{Z) 24z 3 ' 2 36£ 



(6.04) 



The £-map of a suitable ^-progressive path linking C\ with c-i is shown in figure 4. As the radius 
R of the semicircular part tends to infinity, the variation of F along the whole path vanishes. There- 
fore from (5.17) ei(c 2 ) —0. Similarly e 2 (c 3 ) and € 3 (ci) both vanish. Accordingly (5.16) reduces to 

wi(z) + e- 2nil3 w 2 {z) + e 27dl3 w 3 {z) = 0. (6.05) 

By comparing (5.09) with the asymptotic forms of the Airy functions [18], we identify 

m (z) = 27T 1 / 2 Ai (z) . w 2 (z) = 27T 1 ' 2 Ai (ze- 2 ^/ 3 ) , w 3 (z) = 27T 1 / 2 Ai (ze 2 ™/ 3 ). (6.06) 

Equation (6.05) is the well-known connection formula for these functions. It is interesting to ob- 
serve that our theory yields it exactly. 

The vanishing of the error terms in this example is essentially a consequence of there being 
no boundary of D intervening between q and Ck in the neighborhood of infinity. The path of figure 
4 generally cannot be used to establish this in other cases. Using different paths, we formulate 
sufficient conditions for the vanishing of €j(ck) in the next section. 

7. Sufficient Conditions for the Vanishing of the Error Terms 

We again use the notation of section 5, and we suppose, without loss of generality, that the 
turning point zo is at the origin. 

THEOREM 4. Let f(z) and h(z) be regular in the sector S: 71 < arg z < 72 for all sufficiently 
large |z|, and 

f(z) ~Az r , h(z) =0(^ 2 >- 1 " 8 ), as \z\ -> 00, (7.01) 

uniformly with respect to arg z in S. Here A is a real or complex non-zero constant, and r, 8 real 
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constants such that 5 r > — 2, 8 > 0. ,4/so, within S, Ze£ «Sfi, if 2 6e outward-drawn rays from the 
origin which lie in adjacent principal subdomains Di, D 2 , respectively. Then provided that Sfu &2 
are not parallel to the asymptotes of the neighboring boundaries of\}\, D 2 , they are strictly ^-pro- 
gressive paths in the neighborhood of infinity, and if Ci, c 2 denote the points at infinity on them, 
€i(c 2 ) = 0. 

To prove this result, we first observe that 

g(z) = {f(z)}-^~A-^z-'-l* (7.02) 

as \z\ — * °° in S. Hence ([10], section 1.36) 

g"(z) ~ 5 (^+ l) /M^-WO"* (r * 0) (7.03) 

as |z| — * °° in any sector S': yj < arg z < y 2 which is properly interior to S. When r=0, the result 
corresponding to (7.03) is g"(z) = o(z~ 2 ). Therefore, whether or not ris zero, we have 

g(z)g"(z) = 0(z-*- 1 ) (7.04) 

as |z| — ■> °° in S', uniformly with respect to arg z, where 

s — 2 r + 1 , 
and is positive, by hypothesis. From (4.07), (7.01), and (7.04) we derive 

F'(z)=0(z- l -h), (7.05) 

where 8i = min(o\ s). Hence the integral in (5.17) converges. The essential problem is to show 
how to construct a ^-progressive path connecting C\ with c 2 in the neighborhood of infinity. 
From (5.01), we have 

A xl2 z s 
£ = + o(z«) (zGS). (7.06) 

5 

Therefore the principal curves approach infinity in directions parallel to the rays 

arg z= {(n+ 2 )tt— 2 a} Is, (7.07) 

where a = arg A, and ra = 0,±l,±2, .... We call these the principal directions. Consider 
the behavior of Re £ along the ray with parametric equation z = re iy , y being real and fixed, and 
< r < 00. We have 

dildT={f(z)} 1 fHdz/dr)=A 1 l 2 T r f 2 e is y{l+o(l)} (7.08) 

as r— » °°. Hence 

d (Re g)ldr= |^| 1/2 t'-/ 2 {cos (sy-f hot) + o(l)}. (7.09) 

The constant cos {sy+ hot) vanishes only when y has one of the values given by the right of (7.07). 
Therefore any outward ray in S not in a principal direction is a strictly ^-progressive path for all 
sufficiently large \z\. In particular this includes ifi and if 2 . 



Thus the differential equation (4.01) has an irregular singularity at infinity. 
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FIGURE 5. z-plane. FIGURE 6. Approximate path in g-plane: \d\ = \A\ ll2 b/s. 

Next, consider the path with parametric equation 



z={(b + iryy} 118 ' ( - oo < T < oo), 
where y and b are real constants, b being large and positive. 



(7.10) 



As indicated in figure 5 this path begins at infinity, sweeps back to its nearest approach to the 
origin at z— b l/s e ly ^ s and then returns to infinity, the angle between its asymptotes being tt/s. When 
b is large, we have from (7.06) 

£ ~ 5 -i|^|i/2(i + i T ) e f{7+(«/2)} ? 

showing that the £-map of (7.10) is approximately the straight line indicated in figure 6. Cor- 
responding to (7.09), we find that 



1 (Re £ ) = _J^! J sin (y+1 ^ + o(l)} , 



dr 



(7.11) 



as 6— »°°, uniformly with respect to r. Therefore the curve {7.10) is a strictly £ -progressive path 
for all sufficiently large b, except when its asymptotes are in principal directions. 

Now consider c\ and C2. By taking y[ and y 2 sufficiently close to y\ and y?, respectively, we can 
ensure that S' contains if i and if 2. From (7.06) and the hypotheses on £t? \ and J^2, it follows that 



|arg ci~ arg c\\ <2ttJs. 



(7.12) 



Figure 7 indicates how C\ and c>> can be connected by a strictly ^-progressive path in S', consisting 
of if i, if 2 and either one or two arcs of the type (7.10), depending whether or not |arg c^ — arg C\\ 
< tt/s. The corresponding path in the ^-plane is approximately that indicated in figure 8. On 
letting the ^-parameters of the arcs tend to infinity, we see immediately from (7.05) that the varia- 
tions of F along if i and 1£ 2 tend to zero. For the contribution to the variation from the arcs (7.10), 
we observe that the corresponding values of f\z~ 1 ~ 6l dz\ are bounded by 



1 



dr 



|fc + iT |l + <8,/s) 



and accordingly vanish as b—> °°. The proof of Theorem 4 is now complete. 

An alternative set of sufficient conditions for the vanishing of the error terms is provided by 
the next theorem, which can be proved in a similar way. In this result m, A denote real or complex 
nonzero constants, and 8 is a positive constant. 
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FIGURE 7. ^-progressive path in z-plane. 




FIGURE 8. Approximate path in g -plane. 



THEOREM 5. Let f(z) and h(z) be regular in the strip T: 71 < lm(mz) < jifor all sufficiently 
large positive values of Re {mz), and 



f{z) - Ae» 



h(z) = 0(e mz( *-W 2 ), as Re (mz) -+ + °o, 



(7.13) 



uniformly with respect to Im (mz) in T. Also, within T, let ifi, if 2 be straight lines which are 
parallel to the boundaries of T, and lie in adjacent principal subdomains Di, D 2 , respectively. 
Then provided that ifi, if 2 do not coincide with the asymptotes of the boundaries 0/D1, D 2 , they 
are strictly ^-progressive paths in the neighborhood of infinity, and ifci, c-i are the points at infinity 
on them, €1(02) = 0. 

8, The Gans-Jeffreys Connection Formulas (i) 

The first of these two formulas may be regarded as a special ease of the fundamental connec- 
tion formula of section 5. We again consider the differential equation (4.01) in which f{z) and 
h(z) are regular in an unbounded, simply connected domain D and f(z) has just one zero in D, a 
simple zero at 2 = 0. This time we suppose that D includes the whole of the real axis, and that 
f(z) and h(z) are real when z is real. Without loss of generality, we take f{z)/z to be positive on the 
real axis. 




Figure 9. Typical principal subdomains D, for the 
Gans-feffreys formulas. 



283 



One of the three principal curves, say ^i, which emerges from the origin is the negative real 
axis. The other two emerge at directions ±7r/3 and are symmetric with respect to the real axis. As 
in section 5 we denote by Di, D 2 , D 3 the three subdomains bounded by the principal curves which 
have the origin on their boundary. 

We take C\ to be the point at infinity on the positive real axis, and c 2 , c 3 to be points at infinity 
on conjugate ^-progressive curves i^2, ¥z- These curves must lie within D2, D 3 , respectively, and 
be such that no boundary of D intervenes between them in the neighborhood of infinity. We 
assume that |Re £(z)| — > °° when z — » Cj (/=1, 2, 3) (the "dominancy condition") and that the varia- 
tion of the function F, defined by (4.07), converges when z—>Cj. The other assumptions we need 

are 

(i) ^-progressive paths can be found in D which link c 2 or c 3 with any point of the real axis, 

other than the origin, and which do not pass through the origin. 

(ii) € 2 (c 3 ), e 2 (— °°) and € 3 (— °°) all vanish. 

In applications, Condition (i) can be tested by considering the £-map. Sufficient sets of 
conditions for (ii) are to be found in section 7. 

The connection formula (5.16) for the solutions of (4.01) takes the form 



W\ 



( Z ) = e irf/3{i + e 3 (ci)}w2Cz) + e- m73 {l + € 2 (ci)}w 3 (z). (8.01) 



To hx the solutions completely, we define g x (z) to be the branch of {f(z)}~ liA which takes its prin- 
cipal value on the positive real axis (compare (5.08)). The branches of the wj(z) are then interpreted 
in real form as follows. 

When z is positive, we have immediately from (5.09) and (5.10) 

w,(z)=f-^ exp (- [V /2 <fe) {l + €i(z)} (z >0), (8.02) 



|€i(z)|^exp {/ ',.«(F)}-1, (8.03) 

the variation being taken along the real axis. 
When z is negative, we find that 



(8.04) 



Substitution of (5.09) in (8.01) accordingly yields 

M ; 1 (z) = {l + € 3 (c 1 )}e^ 4 |/|- 1 / 4 exp (-i f° |/| 1/2 cfe) {l + e 2 U)} + conjugate. (8.05) 

Let 

1 + e 2 (z) = { 1 + tikz)}e-***\ 1 + e,{z) = {1 + fi(z)}e^\ (8.06) 

where /x(z) and 7){z) are real and /jl(z) ^ — 1. Then by elementary considerations 

\ijl(z)\ ^ |e 2 (z)|, \v(z)\ ^sin - 1 |^(z)| ^tt|€ 2 U)|/2, (8.07) 

provided that |e 2 (z)| ^ 1. Substitution in (8.05) yields 

w 1 (z) = 2{l + f4c i )}{l + fx(z)}\f\- 1 l*smlj° |/| 1/2 <fe + | 7r + r)(z)- rfici)} (z<0). (8.08) 
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From (8.07), (5.17) and (5.10), we deduce that 

\l4ci)\, |2iKci)M^|[exp {i n.r.iF)}-ll (8.09) 

and 

|/x(z)|, |2t;(z)/7t| ^exp {•/ '«,r,(F)}-l f (8.10) 

the variations in both cases being evaluated along f -progressive curves. 

The aggregate of (8.02), (8.03), (8.08), (8.09), and (8.10), may be regarded as the full form of 
the first Gans-Jeffreys connection formula [5], [6], [9]. It expresses the subdominant exponential 
solution on one side of the turning point in terms of the oscillatory solutions on the other side. 
Boundary conditions satisfied by W\(z) are clearly 

€i(z) = 0(l) (z^ + oo); fl(z) 9 7)(z) = 0(l) (*-*-00). (8.11) 

We note that the bounds for the error terms /jl(ci), r)(c\) depend solely on the variation of F along 
the f -progressive path linking c\ and c 2 , and can be minimized by proper choice of path. 

The second connection formula concerns the dominant exponential solution. Unlike the 
subdominant solution, this is not unique. Consider first w 2 (z) . For negative z, we have from (8.04) 
and (5.09) 

w 2 (z) = e-^ 12 \f\- 1 ! 4 exp (- i P \f\ l l*dz\ {1 + e 2 (z)} (z < 0). (8.12) 

On the other side of the turning point 

w 2 (z) = e^y-^ exp ( fV 2 ^) i ] +€ 2 (z)} (z > 0). (8.13) 

Instead of m(z), however, we prefer the real solution w 4 {z) =Re{e- 7ril "w 2 (z)}. Then with the aid 
of (8.06), we derive 

w 4 (z) = [1 + Re{e 2 (z) }]f-i' 4 exp ( |V 1/2 ^) Ifi > 0), (8.14) 

and 

w 4 (z)={\ + fjL(z)}\JV 4 cos ||J l/I^A+^w+^z)} (*<0), (8.15) 

where |e2(z)|, |/x(z)| and \2t}{z)Itt\ are all bounded by the right of (8.10). These relations consti- 
tute the second Gans-Jeffreys formula. We note that w 4 (z) satisfies the boundary conditions 

e 2 (z) = € 2 (ci) + 0(l) (z-> + °o) ; f*(z) 9 ri(z) = o(l) (*-»-»); (8.16) 

where |€2(ci)| is bounded by the right of (8.09). 

9. The Gans-Jeffreys Connection Formulas (ii) 

In this section we attack the problem of section 8 with real-variable theory. Our hypotheses 
are 

(i) The function f(z)/z is positive and twice-continuously differentiable throughout the interval 
— °o <z < °°, and h(z) is a continuous real function in the same interval. 
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(ii) As z^±°o, 



p 

Jo 



f' 2 dz 



is unbounded and the variation of the function F (defined by (4.07)) 



converges. 

In these circumstances it is shown in [15], sections 4, 5 that (4.01) has a solution 6 



«>i(z) = (£//> 1/4 {Ai(D+ 6i(z)}, 



^here 



£-{»■ 



11" w»» 



' 2 dt 



2/3 



|e,(z)| ^k^E-^QMiQiexp {X,f, AH)}- 1] (-oo< z <°o), 



and 



ff(z) = 



24|Cl 



3/2 ± J Ji/i. 



d 2 / 1 



""rfz 2 V|/|'/V l/l 



1 1/2 



<&:. 



(9.01) 

(9.02) 
(9.03) 

(9.04) 



The upper or lower signs are taken in (9.02) and (9.04) according as z is positive or negative, and the 
integration limits in (9.04) are determined by the condition H(0) = 0. The auxiliary functions E(Q, 
M(Q and the constant \ x = 1.430 . . . are defined in [15], section 2. 

On substituting for Ai (Q by means of its appropriate asymptotic forms, we obtain 



M ;i(z) = ^7r- 1 / 2 /- 1 / 4 {l + o(l)}exp (- t'f&dz) (z-» + oo) ? 



and 



^here 



m ( z ) = 7r il2\f\- 



1/4 



sin (Jj/l 1/2 <k + \ *) + PU) + o(l) 



(*-*-«), 



p( 2 ) = 7r l/2| £ |l/4 €l(2r ) 9 



(9.05) 



(9.06) 



(9.07) 



and is bounded. 

In order to throw (9.06) fully into "modulus-phase" form, we appeal to Theorem 4 of [1], This 
shows that solutions of (4.01) exist of the form 



l/l-^U+oOWexp \±i j°\f\ m dz\ as z-»-oo. (9.08) 

Therefore we can express 

^W^TT-^l/l-^ll + ^+o^lsin { f° |/| ^dz + 1 7T + 1/2 + o(l)| (z-*-«) f (9.09) 

where *>i, ^ 2 are real constants to be determined. Comparison of (9.06) with (9.09) shows that 

p(z) = (l + j/i) sin {0(z) + 1/2} -sin fl(z) + o(l), (9.10) 



3 Not to be confused with the W\(z) of section { 
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where 

flto-fVl^+Jw. (9.11) 

Hence p(z) has the form 

p(z) = p/, sin {0(z) 4- i*,} + o(l), (9.12) 

where p/. and ^ are real constants (p/, 2? 0) related to V\ and ^2. Clearly 

p,= lirr7 \p(z)\. (9.13) 

Accordingly, from (9.03), (9.07) and the fact that \t\ l l*M(Q-+ ir 1 !* as {-► — «>, we have 

pz. ^ V [exp {A, '/-'_-,, x (//)} - 1 1. (9.14) 

In (9.10) and (9.12) let z~^ — & through a sequence of values for which 0(z) is an integral multiple 
of 77 (compare hypothesis (ii) ahove). This gives 

(l + i>i) sin Vz = pi d sin r ;? . (9.15) 

Alternatively, by using a sequence of odd integral multiples of 277 we derive 

(1 -h 1^1) cos v-i— 1 = pi. cos v 3 . (9.16) 

Combination of (9.15) and (9.16) yields 

(l-r-^i)e i ^ = l+p / .e i,, 3. (9.17) 

Therefore 

N ^Pu H ^ sin-'p,. ^ 7rp,72, (9.18) 

provided that p L ^ 1 (compare (8.07)). 

The first Cans-Jeffreys formula is represented by the aggregate of (9.05), (9.09), (9.14), and 
(9.18). 

A second solution of (4.01) is [15] 

^t) = (£//) 1/4 {Bi(£)+€2(z)}, (9 j 9) 

where 

|e 2 U)| ^(k 2 lk l )E(QM(Q[exp{\ 1 'r- x , i (H)}-l], (9.20) 

and \> = 1.315 .... Substitution of the asymptotic forms for Bi(£) yields 

w ,U) = 77-" 2 /" ,/4 {l+o1z) + o(l)}exp(|^ (2: -> + »), (9.21) 
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and 



Mz) =77-^ |/| -I/* Jcos(|ui 3/2 +| A + o(l)| (z-»-oo), (9.22) 



^here 



atz) = 7r 1 /^ 1 / 4 exp f-| ^ 2 W). (9.23) 

An obvious extension of Theorem 2 of the present paper shows that a(z) tends to a finite limit, 
cr(°°), as z^ + cg. Accordingly 

w 2 (z)=7T-^f-^{l + o- (oo) + o(l)} exp (J"V /2 &) (*-> + «). (9.24) 

From (9.20) and the asymptotic form of M(£), we obtain 

|o<oo)| ^|V5X 2 Xr 1 [exp{Xi^-oo,»(//)}-l]. (9.25) 

And with the aid of analysis similar to that used for throwing (9.06) into the form (9.09), we readily 
show that 

^(^^^-i^jyj-i^ (i + o(i)} cos [ J l/l 1 / 2 rfzH-i 7rH-o(l)| (*-»-«>). (9.26) 

Relations (9.24), (9.25), and (9.26) consititute the second Gans-Jeffreys connection formula. 

Sharper bounds. Somewhat improved bounds for the error terms of this section can be 
achieved by application of Theorem 3 of [15], i and, in the case of the Bi-type solution, of an in- 
equality of the same character as (3.03) above. The results are expressed by 

p L ^ exp {* '-„oGff)+ X4^o,«(fl)} - 1 ' ( 9 - 2? ) 

and 

|aW|^exp{^^-,,^ < 9 - 28 ) 

where X 4 = 0.521 ... and X 5 = 0.520 .... 



10. Summary 

The primary purpose of this paper has been to investigate approximate solutions of the dif- 
ferential equations (1.01) and, more generally, (4.01) in a complex domain containing no singu- 
larities of the functions j[z) and h(z), and a simple zero of j{z). Some preliminary results concerning 
the Liouville-Green approximation are given in sections 2 through 4. The connection formula for 
the three fundamental approximations of this kind is then established in section 5, complete with 
strict error bounds, by considering limiting behavior at the point at infinity. An illustrative 
application to Airy's differential equation is made in section 6, and some general circumstances 
in which the error terms vanish are investigated in section 7. 

An important special case concerns real variables: the problem here is to connect the exponen- 
tial-type solutions on one side of the turning point with the oscillatory solutions on the other. 
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This is solved in section 8 by use of the general formula of section 5. An alternative solution is 
constructed in section 9 solely from real-variable theory, by re-expansion of known uniform Airy- 
function approximations. Naturally, different hypotheses on J(z) and h(z) are made in the two 
methods, but the principal difference lies in the forms of the error bounds. In the first method 
the bounds depend on the variation of a certain function along an infinite path of specified type in 
the complex plane. In the second method the bounds depend on the variation of a more compli- 
cated function along the real axis. 

The error bounds given in this paper have computational applications. As shown in the paper 
which follows [19], they may also be used to provide theoretical insight into the accuracy of the con- 
nection formulas in various circumstances, including cases when the differential equation contains 
a large parameter. Hitherto, in certain applications it has often been difficult to assess the validity 
of assumptions which have been made. 

11. Addendum 

The problems studied in the present paper and its sequel [19] have also been considered in 
a recent monograph by Froman and Froman [20]. These authors have the same principal objec- 
tive, that is, the determination of strict error bounds for the coefficients in the connection formulas. 
The method they use resembles the first method of the present paper in that it is a rigorous formula- 
tion of the procedure of Zwaan. The difference in the approaches is as follows. 

We constructed (in sec. 5) the connection formula for the three fundamental Liouville-Green 
solutions associated with a simple turning point of the differential equation (4.01) by examining 
the asymptotic behavior of these solutions in the neighborhood of the point at infinity. Froman 
and Froman, following Kemble [3], express the general solution of (4.01) in the form 

w{z) =a l (z)f- 1 ' 4 exp ( f fl* dz) + a*(z) f" 1 ' 4 exp (- J/ 1 / 2 dz), 

and derive a pair of first-order simultaneous differential equations for the functions a\{z) and a-ziz). 
These equations are then used to determine the changes in these functions as z passes from one 
principal subdomain to another. Error bounds for the changes are expressed in terms of contour 
integrals which are equivalent to the variations of our function F(z). 

Generally the path for the contour integrals consists of the real z-axis with indentations by- 
passing the turning points. Froman and Froman show that the appropriate form of monotonicity 
condition is fulfilled when the indentations are semicircles of sufficiently small radius ([20], p. 
39). Owing to the singularities of F{z) at the turning points however, the error bounds tend to 
infinity as the indentations shrink to zero. In consequence Froman and Froman assume that the 
indentations can actually be taken sufficiently large to yield small error bounds, without violating 
the monotonicity condition. In the present paper we allow greater flexibility in the choice of paths, 
with a view to achieving minimum bounds in applications. 

Other comparisons with the first method of the present paper are as follows: 

(i) There is greater emphasis in [20] on the case in which the coefficients f(z) and h(z) in the 
differential equation (4.01) are real on the real axis. 

(ii) Approximations and expansions are frequently used in [20] to simplify the error bounds. 
There is greater emphasis on strict bounds in the present paper and [19]. 

(iii) In the special limiting case in which z is at infinity the error bounds can be sharpened 
by a factor of 2 (see Theorem 2 of sec. 3 above). This refinement is freely used in the present paper 
and [19], but is not in evidence in [20], 
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